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Lattice calculations could boost our understanding of Deep Inelastic Scattering by evaluating mo- 
ments of the Nucleon Structure Functions. To this end we study the product of electromagnetic 
currents between quark states. The Operator Product Expansion (OPE) decomposes it into ma- 
trix elements of local operators (depending on the quark momenta) and Wilson coefficients (as 
functions of the larger photon momenta). For consistency with the matrix elements, we evaluate 
a set of Wilson coefficients non-perturbatively, based on propagators for numerous momentum 
sources, on a 24 3 x 48 lattice. The use of overlap quarks suppresses unwanted operator mixing 
and lattice artifacts. Results for the leading Wilson coefficients are extracted by means of Singular 
Value Decomposition. 



The XXVI International Symposium on Lattice Field Theory 
July 14-19, 2008 
Williamsburg, Virginia, USA 



* Speaker. 



© Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike Licence. 



http://pos.sissa.it/ 



Nucleon structure in terms of OPE with non-perturbative Wilson coefficients 



W. Bietenholz 



1. Motivation 

The computation of moments of the nucleon structure functions is a fascinating challenge: it 
is known to be difficult, but it is a point where lattice results could contribute much to the under- 
standing and interpretation of phenomenological data from Deep Inelastic Scattering. 

The continuum formulation of hadron structure functions is plagued by renormalon ambigui- 
ties, i.e. power-like IR contributions, see e.g. Refs. Here we refer to the lattice regularisation, 
where a general moment of the nucleon structure function can be expanded as 

Jt(q 2 ) =c^ 2 \aq)A2(a)+c ( - 4 \aq)\A 4 (a) + ... { higher twists } . (1.1) 

q 1 

It depends on the transfer momentum q, while a is the lattice spacing, c^ , c^ ... are Wilson coef- 
ficients (where the superscript is the twist), and A2, A4. .. are reduced matrix elements ("reduced" 
in the sense that the Lorentz structure is factored out). 

For the evaluation of A2 there is an established procedure, which employs the ratio between 
two-point and three-point correlation functions [Q]. On the other hand, the Wilson coefficients 
have usually been evaluated in continuum perturbation theory. However, we need a cancellation 
of singularities in the terms c^ , A4 oc \/a 2 , which are again a facet of the renormalon problem. 
This requires a strictly consistent treatment [f||. Therefore we evaluate the Wilson coefficients 
non-perturbatively as well. This method is particularly suitable for disentangling higher twist con- 
tributions. The use of overlap quarks suppresses undesired operator mixings. 

In this report we present new numerical results for c^ 2 \ For earlier results with Wilson 
fermions we refer to Ref. [^j. In a previous study with overlap quarks on a 16 3 x 32 lattice [||] 
some problems persisted, which motivated us to enlarge the lattice to the size 24 3 x 48. 

2. Operator Product Expansion on the lattice 

To be explicit, we apply the Operator Product Expansion (OPE) to a product of electromag- 
netic currents between quark states \\ff(p)), 

W^ v (p,q) = (Y(p)\Jn(q)4(q)\Vf(j>)) 

m 

The lower line is a decomposition in terms of local operators tf( m > , which characterise the nucleon 
structure. The index i = 1 ... 16 specifies the Dirac structure, pL\ ...\l n indicate the momentum 
components involved, and m distinguishes operators with the same symmetry. Note that the corre- 
sponding Wilson coefficients C^ m > depend solely on the transfer momentum q. 

A truncation of this expansion at some low operator dimension requires the following scale 
separation, 

p 2 «<? 2 «(^) 2 . (2.2) 

So far we have been using periodic boundary conditions for the gauge field. Thus a sizable lattice is 
required to have a set of small momenta p 2 available. If this scale separation is realised, it justifies 
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a truncation of the OPE (2.1 ), which cuts off operators with high derivatives, i.e. high powers of p. 
Here we consider quark bilinears with up to 3 derivatives, 

YTiY, yTiD^y , yTiD^D^y , yFiD^D^D^y , (2.3) 

where T, runs over a basis of the Clifford algebra. This still amounts to an apparently frightening 
set of 16 ■J^ =0 4 d = 1360 operators. However, we choose the isotropic transfer momenta 
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(| gB |~4.1GeV), 
(\q b \ ~5.5 GeV) . 



(2.4) 



In this case the symmetry reduces the set of operators to only 67 equivalence classes We denote 
them as C\ . ..Cei, which refer to an ascending number of derivatives: 

C\ : no derivative, vanishes in the chiral limit, multiplies yly 
C%...C(, : one derivative, Bjorken scaling ^ l/q 2 
Cj . . .Ci6 : two derivatives, vanish in the chiral limit 
Cn . . .Ctf : three derivatives, Bjorken scaling « 1/ (q 2 ) 2 ■ 

The coefficients of terms with an even number of derivatives vanish at quark mass m = due 
to chiral symmetry. In the Bjorken limit of large q 1 the coefficients of terms with one (three) 
derivative(s) are expected to scale as C m l/q 2 (C m °< l/(q 2 ) 2 ). 

To separate the scales even better, we have now implemented twisted boundary conditions: 
thus very small p 2 become accessible, which enables us to use q 2 further below the momentum 
cutoff squared [§]. This allows us to consider for instance q = ^(1,1,1,1). 

3. Results for the Wilson coefficients 

In our numerical study we analysed configurations that were generated in the quenched ap- 
proximation with the Luscher-Weisz gauge action in a volume V = 24 3 x 48, with a physical lattice 
spacing a ~ 0.075 fm. To provide finite OPE matrix elements, as they occur in eq. (|2~I|), we fixed 
the lattice Landau gauge; this is also favourable to reduce the statistical noise. 

For the valence quarks we used overlap fermions with the parameter p = 1 .4 (negative mass 
of the Wilson kernel) at a bare quark mass of 0.028 in lattice units, which corresponds to about 
73 MeV. The application of chiral lattice fermions suppresses undesired operator mixings, as well 
as 0(a) lattice artifacts. 

We denote the number of /j-sources as M, and the number of Wilson coefficients to be con- 
sidered as C\ . . .Cn, with N = 61 in our case. So far we evaluated W^ v off-shell for M = 25 
/^-momentum sources at transfer momentum q a (on different configurations), and for M = 10 p- 
momentum sources at qt,. In each case this yields a system constraining the Wilson coefficients, 
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(3.1) 
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where the elements and & k are 4 x 4 matrices. 
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Figure 1: The convergence of the Wilson coefficients as the number of the n most powerful constraints is 
increased. We show the residues as a function of n for the transfer momentum q a with 25 p-sources (upper 
plot), and for q/, with 10 ^-sources (lower plot). In both cases we see a convincing convergence before n 
reaches the deterministic number 67, which confirms that the extracted Wilson coefficients are trustworthy. 

In both cases the system is over-determined since 16M > 7Y. We apply Singular Value Decom- 
position as an established method to analyse such systems 0. Roughly speaking, this methods 
selects the n < N conditions with the "maximal impact" on the solution C\,... ,Cjy. A rapid con- 
vergence for increasing n approves a reliable result {i.e. the remaining conditions are negligible). 
With 12 singular values (analogues of eigenvalues in rectangular matrices [^]) we do observe this 
feature. As an example Figure [j] shows the convergence in n for the Wilson coefficients under con- 
sideration, which saturates around n w 50 for q a and for qb- In the latter case, the total number of 
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constraints is lower, hence it is easier to satisfy them to a good precision. Therefore the residues 1 
are smaller, but the result is less reliable. Note also that q 2 = \(n/a) 2 is more promising than 
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^(n/a) 2 in view of the required scale separation (2.2). 
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Figure 2: The Wilson coefficients evaluated at momentum ttansfer q u (upper plot) and qt, (lower plot). The 
results are compared to the corresponding tree level values (with the same overlap Dirac operator). We 
observe the same pattern, but a significant non-perturbative correction. In the chiral limit C\ and Cj . . .Ci6 
vanish. The small values that we obtain (with light quarks) for these coefficients provide a consistent picture. 



The residue is the norm of the difference between the two sides of eq. (3.1 ) in the given approximation 
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There are a number of discrete rotations and reflections which leave our isotropic ^-momenta 
(2.4) invariant. For example this allows us to exchange W33 with W44. As a consequence, specific 
pairs of Wilson coefficients belonging to different operators have to coincide; for instance the co- 
efficient of ^ 33 in W33 is equal to the coefficient of ^44 in W44. We implemented this property in 
combined fits, which have of course somewhat larger residues, see Figure [[]. 



We now proceed to the actual results for the Wilson coefficients obtained from W33 and W44. 
The best results emerge from combined fits. They are shown in Figure ||, which also displays the 
corresponding tree level values for comparison. 

In the limit of zero quark mass the coefficients C\ and Ct ... C16 vanish due to chiral symmetry, 
as we anticipated in Section 2. Since we are dealing with light quarks represented by chiral lattice 
fermions, we are fairly close to chirality. Therefore it is a stringent consistency test that we do 
obtain particularly small values for these coefficients. That property had not been observed with 
Wilson fermions [Q]. 

Finally we also test the Bjorken scaling property, which we mentioned in Section 2. We show 
again the C„, obtained at q a . We compare them to the Wilson coefficients evaluated at qb, where 
C2 . . -C(, are enhanced with a factor q\lq a = 16/9, and C17 . . .C(,j are amplified with the square 
of this factor. The Bjorken re-scaled results are confronted in Figure |[ and we see an impressive 
similarity, in particular for cases with relatively large \C m \. 
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Figure 3: The Wilson coefficients evaluated at momentum transfer q a and qb, where in the latter case 
the coefficients of terms with one or three derivatives are amplified according to the expected scaling in 
the Bjorken limit. In fact the C m are now in the same range, and we observe a very satisfactory level of 
agreement. 
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4. Conclusions and outlook 

Our method to evaluate the Wilson coefficients non-perturbatively works successfully. We ar- 
rived at conclusive non-perturbative results for a considerable set of Wilson coefficients. We are 
continuing to evaluate constraints for additional q and p momenta to further tighten the control 
over possible artifacts. The use of twisted boundary conditions enables us now to achieve an even 
better scale separation. 

In a final step, the nucleon structure function ^# is obtained by Nachtmann integration [||] 
over Way. For instance, for the second moment this integral takes the form 



•^2(0) = / da <l n li Wn V (q)-~8n V Wpp 

f l r ? 1 

-> / dx F 2 (x,q ) + -F L (x,q' 
Jo L 6 



(4.1) 



where the Bjorken limit is taken in the lower line. The projection vector has length n 2 = 1. For a 
different projection one obtains J dx [F2 — |fzj instead. Thus the combination of different projec- 
tions determines the longitudinal structure function Fl = F 2 — 2xF\. 

With completed data, we will obtain a fully non-perturbative and consistent evaluation of the 
Nucleon Structure Functions [p. 
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